Andreev reflection in superconducting QCD 



Mariusz Sadzikowski 1,2 and Motoi Tachibana 3 
1) Center of Theoretical Physics, Massachusetts Institute of Technology, Cambridge, MA 02139, USA 
2) Institute of Nuclear Physics, Radzikowskiego 152, 31-342 Krakow, Poland 
3) Theoretical Physics Laboratory, RIKEN, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan 

MIT-CTP-3289 



< 



X 



In this paper we discuss the phenomenon of the Andreev reflection of quarks at the interface 
between the 2SC and the Color-Flavor-Locked (CFL) superconductors appeared in QCD at asymp- 
totically high densities. We also give the general introduction to the Andreev reflection in the 
condensed matter systems as well as the review of this subject in high density QCD. 



I. INTRODUCTION 
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Quantum Chromodynamics is a non-abelian quantum field theory describing interactions between quarks and 
gluons. It describes the effect of confinement at low energy [Q and asymptotic freedom at high energy Q compared to 
the QCD scale Aqcd ~ 200 MeV. The only known systematic approach to low energy QCD are lattice calculations. 
The physics of the high energy QCD can be controled by the perturbation analysis in strong coupling constant. There 
are three different types of phenomena of strong interactions which one can hope to describe in the perturbative 
\ regime: high energy scattering processes, high temperature and high baryon density systems. The energy scale for 
these phonomena are set by: the value of the four-momentum transfer •*/ Q 2 exchanged in the process, the value of 
f — , temperature T and the value of the Fermi energy Ef (or quark chemical potential fx) respectively. The high energy 
scattering processes were already well tested in deep inelastic scattering of leptons on protons or electron-positon 
annihilation (DESY, CERN) . The samples of the matter at high temperature (T ~ 200 — 300 MeV) and low baryon 
density one can hope to achieve in high energy ions colliders (RHIC). However lattice calculations show that the 
expected temperatures are too low by the factor of 3-4 for the application of the perturbative QCD. Then one has to 
rely on the non-perturbative analysis: lattice calculations, universal features of phase transitions or models. FinaUy 
■»«. ' the high baryon density quark matter at low temperatures is expected to exist in the cores of the compact stars |p| . 
'■^ Unfortunately, here also the densities are not high enough for direct perturbative calculations (by many orders of 
I 1, magnitude). Then the only possible approches are based on universal features of the phase transitions and models. 
Q_i| The lattice approach fails at non-zero density physics because of the sign problem. 

In this paper we focus on the high baryon density phases at T = 0. Since long ago one expected supeconducting 
phenomena in this kinematical region Q. However only recently there have been a new interest in the subject because 
new features were found in the theory |^ M . If the energy scale is set by the quark chemical potential fx one can 
expect that at high enough densities there are a gas of free quarks and gluons. However the interaction between quarks 
mediated by the one-gluon exchange is attractive in the color 3 channel. This leads to the well known phenomenon of 
■ the Cooper instability of the Fermi sea and finally results in the creation of a new vacuum - the condensate of Cooper 
pairs. This is the version of the BCS theory of superconductivity || applied to the interacting quark matter. The 
picture is generally correct, however somehow oversimplified. There are subtle but important differences between BCS 
theory and superconductivity in QCD which follow from the non-abelian character of the interation in quark-quark 
scattering (eg. 0). The review of the whole subject can be found in j|. Despite of some differences the high density 
QCD shares a lot of features with condensed matter systems. Indeed at high density the relativistic quantum field 
theory effects are suppressed and one can apply just non-relativistic field theory or relativistic quantm mechanics to 
study interesting phenomena. 

Our subject of interest are scattering processes at the interfaces between different types of QCD phases. These 



are the analogs of the Andreev reflection |10| in condensed matter systems. This reflection apperas at the junction 
between conductors and superconductors. In the case of dense QCD the role of conductor is played by the free Fermi 
gas of quarks and superconductor is 2SC or CFL phase of QCD. In this paper we consider the interface between 2SC 
and CFL phase. We shall show that this interface is the most general in a sense that it already contains the other 
cases: free quarks/2SC fill and more: 2SC/2SC and one only peculiar to 2SC/CFL. The case of free quarks/CFL 
interface has to be considered independently. 

The interest in the dynamics of the matter flow through the phase boundaries is of importance because such 
interfaces can be present in the protoneutron stars. Indeed in the final stage of neutron star cooling one expects that 
the star can go through the all possible phase transitions |l3[ including 2SC/CFL case. This last transition is expected 
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to be the first order |T4j thus one expects the existence of the mixed phase: CFL-bubbles in the 2SC medium as well 
as 2SC-bubbles in the CFL medium. In such the conditions the Andreev reflection is frequent phenomena and should 
be taken into account in the calculations of detail processes of the protoneutron star evolution. 

In the sections 2 and 3 we provide the general introduction to superconductivity and Andreev reflection in condensed 
matter systems. Section four contains short description of 2SC and CFL phases, whereas in Section five we describe 
the Andreev reflection in QCD. There are also three Appendices which gives more detailed calculations. 

II. CONDUCTORS AND SUPERCONDUCTORS 

Let us start with the simple effective hamiltonian describing superconducting state of the system at zero tempera- 
ture: 



H = / d A x 



J2 ^.(*>**) - Ef ) Mt,r) + A(r)^\(t,r)^\(t,r) + A^r^^rj^i^r) 



(1) 



The fermi field operator ip a (t,r) (^a(t,r)) describes the annihilation (creation) of the particle of mass m and spin 
a at the point r and the time t. The particles also carry charge however in our simple example we do not consider 
the electromagnetic fields thus this feature is unimportant for our purpose. The quantity Ep defines the Fermi 
surface below which all energy levels are occupied in accordance with Pauli principle. The so called pair potential 
A(r) varies in space and describes the superconducting gap in excitation spectrum of the system - the essence of 
superconductivity. The last 2 terms in the above effective hamiltonian are a mean-field approximation of attractive 
four-fermi point interaction. The pair potential A is related to the anomalous particle-particle correlator: 

A a/3 (t,r) = (V> a (*,f)#(*,r)) = A (t,r)e al3 . (2) 

One can treat this hamiltonian as the description of the interation between particles ip and given classical pair potential 
A(r). We assumed that this interaction is independent of the spin. In that sense the spin is not dynamical degree of 
freedom but only a quantum number that distinguishes between 2 types of particles. 

Let us mention that the hamiltonian (Q) does not conserve particle number^. This results from the fact that the 
pair potential describes the condensate of Cooper pairs which is a source of particles in our approximation. Let us 
also notice that the energy is conserved however not necessarily momentum. The momentum conservation is related 
to the space dependence of the gap parameter A(r). 

The equations of motion that follow from the hamiltonian (| l|) a nd standard anticomutation relations for fermi fields 
take the form of so called Bogoliubov - de-Gennes equations ]fL5| : 



idt Ui^) = { aV) & +Ef (3) 



These operator equations are linear which means that we are free to treat them as the usual wavefunction equations 
(hermitian conjugate changing to complex conjugate). The physical interpretation of the wavefunctions comes from 
the following reasoning. In the case of A = the equations (Q) decouple from each other: 

*Vr (*> »0 = (- ^ - e f) V>T (*, (4) 
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The plane wave solutions: 



( <M*>^ \ = ( f exp(-i£< + iqr) 



(5) 



1 In this situation charge is also not conserved. 
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lead to the dispersion relations: 



E = e q for particles ,., 
E = —e q for holes, ^ ' 

where e q = q 2 /2m — Ep is a measure of the kinetic energy with respect to Fermi level and E describes the energy 
excitation above (below) Fermi energy for particles (holes )Q. Let us mention the interesting detail that for the holes 
the group velocity is opposite to momentum v = —q/m thus the wavefunction ip*(t,r) describes the hole which is 
moving from right to left for r = (0, 0, z)^. The above physical picture can be treated as a simple description of the 
conductors. The energy spectrum contains particles and holes in the vicinity of Fermi surface which can be excited 
easily for the negligible cost of energy. 

Let us now consider equations (|3|) with A(r) = const., which gives us the model of superconductor. The plane- waves 
are still the solutions of (0): 



(7) 



f q e~xp(—iEt + iqr) 
g q exp(—iEt + iqY) 

where momentum dependent constants satisfy the equations: 

- B £-J- „) (;:)=»- 

The selfconsistency condition for the set of homogenous equations gives the dispersion relation: 

E 2 =e 2 + \&\ 2 . (9) 

This is famous gapped energy spectrum of the quasiparticles in the superconductor. The wavefunction describing 
quasiparticles are given by the solutions of (ph: 



J±(l+£/E)ex.p(i6/2) \ ( J 1(1 - (,/E) exp(i£/2) 

= D\ 1—' yy 1 ' ) exp(i<f + r — iEt) + F I VJ± — — 1 ' ' I exp(iq_f- iEt), (10) 

i(l - i/E) expHJ/2) / I J\{1 + i/E) exp(^6/2) 



where ^ = \J E 2 — | A| 2 , 8 is a phase of the gap parameter, D, F some constants and momenta of the excitations are 
given by the relations: 

|j-Sr±{. (ID 

The first (second) wave- function of ([[o|) is particle- like (hole- like) excitation because for E » |A| it reduces to the 
particle (hole) wavefunction. Let us notice that for the excitation with energy E < |A| the momenta are complex and 
thus quasiparticles do not propagate in the medium in such a case. 



III. ANDREEV REFLECTION IN CONDENSED MATTER SYSTEMS 

Let us consider the conductor - superconductor juction (NS). This can be modeled by the equations (0) when the 
pair potential takes appropriate shape. In the simplest case one can consider the plane junction perpendicular to 
the z - axis located in the point z = 0. Then the gap parameter is only a function of z coordinate A(r) = A(z). 
The exact form of the pair potential depends on the details of the junction however some general features can be 
already obtained in the simple model of the step junction A(z) = A6(z) which describes a conductor for z < and 
superconductor with a constant gap A for z > 0. The task is to solve the stationary problem of coupled Schrodinger 
equations (JsJ) with boundary conditions: 



2 The apparent correlation between spin and particle-hole interpretation is accidental. We can have the same set of equations 
with exchanged spin indices. 
3 One can also say that the hole has a negative mass. 
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1. For z — » — oo excitations are particles or holes (conductor) 

2. For z — > +00 excitations are quasiparticles (superconductor) 

3. For z = wavefunctions and their first spatial derivatives are continuous. 

For definitness let us consider that the particle of given energy E incoming from the left hits the junction at z = 0. 
After the collision the particle and/or hole can be reflected back and the qusiparticles can be excited in superconductor. 
Thus for the conductor side (z < 0): 



exp(ikz) + B exp(-ikz) 
Cexp(ipz) 



(12) 



where k — y/2m(Ep + E), p = y/2m(Ep — E) and the constants B and C are related to the probability of the 
reflection for the particle and hole, respectively. For the superconductor side (z > 0) the wave- functions are given by 



formula ( |i0[ ) with f = (0, 0, z). The momenta are given by: q± — ±to + im\£\/w where w — y m(Ep + yj Ep — £ 2 ) 

for E < A and q± — ±y/2m(EF ± £) for E > A. The constants D and F describe the probability of the particle- like 
and hole-like quasiparticles excitations. The continuity conditions at z — lead to the 4 equations for 4 constants 
B, C, D and F. The solutions of these equations are: 

b = ° { e\ ] (13) 



Thus one can conclude that the dominant contribution comes from the hole reflection in conductor and particle- 
like excitation in superconductor. The other contributions are suppresed by the powers of Fermi energy. For the 
interpretation of our result let us consider the transition coefficients. From the equations ([}]) follows the conservation 
of probability current: 

^(IV'tI 2 + IV'il 2 ) + v • ^ (vfvvt - VnWf + ViVVI - ^tvvi) = o. (w) 

This probability current through the NS junction has a form: 

. _f0 for£<|A| 

Jz ~\^F+0(^) for J B>|A| ^ 



where vf — \/2Eplm is Fermi velocity. The result of the vanishing current for energies below the gap is exact. The 
transition and reflection coefficients read: 



Rhole = 1, T quasl = for E < |A| (16) 
Rhole — g ; 7 ; T quasi = =—— for E > |A 



Thus for the energies of the incoming particle below the gap only the hole is reflected. This is easy to understand 
because there is no possibility to transport energy into the superconductor below the gap. Let us notice that in the 
Andreev reflection the charge is not conserved. This follows from the fact that our hamiltonian does not conserve 
particle number. Obviously at the microscopic level the charge is conserved. The interpretation of the process is 
simple: the incoming particle creates the Cooper pair with another particle from the Fermi sea and "dilutes" in the 
condensate. The lack of the second particle is visible as a hole. Additionally the Andreev reflection process conserves 
the energy and momentum (up to 1/Ep corrections). It is also called retro- reflection because the outgoing hole after 
reflection is a time reversal picture of the incoming particle before the reflection. This is approximate picture exact 
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only in the limit of the vanishing particle energy E — 0. Although the probability current vanishes at the NS junction 
the charge current is excessed. Indeed the charge current at the conductor side is double because the particle and hole 
currents add to each other. The charge current is conserved at the interface thus at the superconducting side we also 
have to encounter doubled charge current. This current is carried by the condensate (supercurrent) itself because the 
quasiparticles penetrate only on the depth of the order vf/|£I- 

For the energy of the incoming particle above the gap the coefficients Rhoie and T quas i give us the probabilities of 
the reflection of the hole or the transmission of the quasiparticle into superconductor. For E » |A| the coefficients 
Rhoie — > 0, T quas i — > 1 as it should be. 

These predictions were checked experimentally many times, e.g., by the direct measurement of the reflected holes 
fL6|| . The Andreev reflection can influence transport processes at the NS junctions which were already shown by 
Andrccv jlC| ]. Other effects like excess current, charge imbalance and supercurrent conversion were also considered 
(e.g., @). 

IV. SUPERCONDUCTIVITY IN QCD 

The perturbative approach to high density QCD is well established at much higher densities than one can expect 
in the cores of the compact stars. Thus the use of models are inevitable in the description of interesting physics. This 
is the approach we chose in this paper. The one gluon exchange suggests that the diquark condensate is created in 
color and flavor anti-symmetric scalar channel. Similar feature one can also find using instanton-based models. The 
effective hamiltonian describing quark interaction with the pseudoscalar condensate at the high baryon density can 
be written in the general form (for the review see [[!): 



H = d 



j3„ 



^t H5 . v _ M) ^ + y Arwfc 7 ^) + ii. c 



a,i a,b,i.j 



(17) 



where tp a are Dirac bispinors, a, b color indices, i,j flavor indices, C the charge conjugation matrix and \i is the 
quark chemical potential. At zero temperature the condensate A A is a vacuum expectation value of the 2-point field 
correlator: 

WaC^l) = A^, (18) 

The formula ( |l7| ) is a relativistic generalization of the hamiltonian (Q) with the more complicated gap structure. 
In the case of two flavors, the gap matrix ( |l8| ) becomes: 

A« = Ae^e ab3 (19) 

where the third direction in color space was chosen arbitrarily. This combination breaks gauge symmetry SU(3) C — > 
SU(2) C whereas chiral symmetry SU(2)l x SU(2)r remains untouched. This is the 2SC phase. The Cooper pairs are 
created between the quarks of different two colors and flavors and third color quarks are unpaired. The lowest energy 
excitations are of course unpaird quarks. The quasiparticles are separeted by the gap A from the vacuum. The value 
of the gap depends on the model and is usualy in the range 50 — 150 MeV. 
For the case of three flavors the gap parameter takes the formQ: 

A^ = Ae» k e kab . (20) 

This is color-flavor locked (CFL) phase. The gauge symmetry and global chiral symmetry are broken according to 
the scheme SU(3) C x SU(3)l X SU(3)r — > SU(3)l+r+ c - The baryon U(l) symmetry is also broken. All quarks are 
gapped in Cooper pairs with different colors and flavors. The lowest excitations are nine Nambu-Goldstone bosons 
related to the spontaneous symmetry breaking of global symmetries. There is also additional Nambu-Goldstone boson 
connected to the breaking of restored axial symmetry [/^(l). The value of the gap is model dependent of the order 
of the 2SC gap but usualy slightly smaller. 
The gap matrix can be written in general as: 



4 There is also a small admixture of the condensate in color symmetric 6 channel but we neglect it in our considerations. 
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/ 

A ud 



A ab 

13 



A u d A us 

A ds 
A ds 





-A u d 



-A u d 














-A„ 



-A ds 




(21) 



in the basis 



where: 



i^red 1 d green 1 $blue 1 ^red 1 ^ green 1 $red 1 ^blue 1 $ green 1 ^blue)^ 



A rfs 



0, A ud = A for 2SC 
A ud = A, for CFL 



(22) 



(23) 
(24) 



From the gap matrix structre one can easily recognize the (u re d,d green ), (u 



L green> ^re 



d) pairings in 2SC phase. The 

quasiparticle excitations are two doublet representations of the remaining unbroken symmetry SU(2)l X SU(2) C 
of given chirality L. The second similar set exists for oposite chirality. The unpaird quarks are singlets un- 
der color rotations SU (2) c and doublets under chiral symmetries. All of the excitations have the same gap 
A. The detail wave- functions of quasiparticle are given in Appendix B. For the CFL phase we have pairings: 



(^retZj dgreen) j ($redi ^blue) > 



green ) 



d-biue) similar to the 2SC and the combination of (u re d, d 



green 7 $blue 



which is essen- 

tialy new for the CFL phase. After diagonalization of the matrix we can find 8 excitations with the gap A and one 
excitation with the gap 2A. These are octet and singlet representations of the unbroken SU(3)l+r+ c symmetry. The 
detail calculations are given in Appendix C. 



V. ANDREEV REFLECTION IN SUPERCONDUCTING QCD 

In this section, let us consider the Andreev reflection at the interface consisting of two different superconductors 
in QCD, that is, the 2SC/CFL interface since we find it is the most general case among the possible interfaces. As 
is obvious from the expression of the gap matrix (21) with the conditions (23) and (24), there are three kinds of 
possibilities of the quark scatterings at this interface. 

• (sredjUbiue) and/or (s green , duue)'- in this case, the Andreev reflection is similar to QGP/2SC because strange 
quark as well as blue up and down quarks are unpaired in the 2SC phase. This interface has already been 
studied Q. 

• (d red , u green ): in this case, the reflection is 2SC/2SC - like because the gap of each phase is, in general, different. 

Sbiue)- in this case, the reflection is a mixture of (QGP-2SC)/CFL interface. This happens because 
in the 2SC phase the u re d, d green quarks are paired whereas strange quarks remain free. 

Since the last two cases are specific to the 2SC/CFL interface and have not been studied so far, let us restrict our 
consideration into those cases. The comparison to other possibilities would be given at the end of this section. 

Let us start with the 2SC/2SC case. Physical setup we are interested in is that the interface is placed at z = in 
space. For z < 0, we have 2SC phase and for z > 0, CFL phase and we provide the boundary condition that matches 
the wavef unctions at z — 0. 

For z < 0, the wavefunction takes the form as (below we put m = for simplicity) 



E-X h djT 
2E n iL 



G 



while for z > 0, 



*>(*) = C 



e 2 



2B ^T-R 



2£ n iL 



D 



2£ "lL 



(26) 



where fci = + A, fc 2 = — /i + A,pi = p + £ and p 2 = — A* • £• i \JE 2 — A 2 and A = \J E 2 — A 2 . A is the gap in 
the 2SC phase. The angles b s and S c are phases of the gap parameters in 2SC and CFL superconductors respectively. 
The bispinors ipn, are desribed in Appendix A. 

By matching these wavefunctions at z = 0, we find the results as follows; 



A = 0, 
B = 



C 
D 





0(E- 


-xy a -^(E + 0(E-\) 






-A)-V(^-0(^-A)e- 






2Ae i f 






-A)- V(^-0(^-A)e iQ 



0(1/ M ), 
0(1///), 



0, 



(27) 



where a = S s — S c is a phase difference of the gaps crossing the interface. Note here that the coefficients A and D 
exactly vanish in the massless limit. First of all, the transition and reflection coefficients, when E > A, A, are given 
by the formulae: 



T Lg\* ^ 1 

A 1 + + A) 1 + r 2 -2rcosa' 

R = l— T 



where 



'(E-\)(E-0 



whereas for A<£'<A,T = and of course R 
coincides at both sides of the interface: 



/ (E + \){E + £) 

1. The conserved current is defined as j 
for A > E > A 



(28) 



(29) 



ip < dip < and 



J 



2 ^E (B+C)(iS+A)(l+r ;i -2rcosa) fol ' E > A, A 



(30) 



The coefficients and currents depend on the phase difference a of the gap parameters. We suggest that this combination 
is a gauge invariant quantity. Indeed in the case of U(l) superconductors a variable is responsible for many interesting 
physical phenomena like, for example, the Josephson effect. 

Let us move to the case (QGP - 2SC)/CFL interface (the mixed case). In this case, the basis of the wavefunction 
must be 



For z > 0, we can use the equation 



■ <5 S 

e 1 2 



E+X u 
IE r]R 







( r red \ 

ib d 

Y green 

^blue 

■ red 
JjdfT 
Y green 

|) from Appendix C. In the 2SC side the wavefunction takes the form: 



(31) 



Ak\ z 



A 



IE r]R 







E-X ^dfT 
IE n lL 





-ik\z 



B 



e 1 2 



2E r]R \ 







2E ,l lL 





7 



+ c 



V 







E+X d 
2E r]R 





E-X h ujT 
2E n lL 



— ikiz 



D 



J 



■ 8a 

e 1 2 



V 



o 



2E r-[R 





E+X h «t T 
2E n lL 






/ 



-ik,2Z 



( \ 









V o J 



— ik%z 



G 



( \ 






\ /l sfT / 



-ik^z 



(32) 



The first wavefunction describes the incoming quasiparticle, second and third the reflected quasiparticles of the same 
kind. The next two wavefunctions are quasiparticles of the second kind whereas the last two describe the posssibility 
of the reflection of the unpaired strange quark. The momenta of the strange quark and hole are = [i + E and 
&4 = —fx + E, the other were already given previously. On the other hand, for z > 0, we have the wavefunction given 
by the formula ( |54| ) in Appendix C. In the massless limit we consider here the momenta in the wave function ( |54| ) 
are given by q\ = fi + £, q2 = —li + £, pi = fj, + £ and P2 = —(J, + (■ Solving the boundary condition which connects 
^<(z) with ' l S > (z), we find the following results (the phase exp(i8 s /2) was absorbed into the redefinition of constants 
G,H,K and iV): 



B = -l 



D = - 



2e la {e la lx{x - 2z) + x + z)X 
(1 - e ia lx)(3 + e ia x(x - 2z))(E + A) 

2e ia {2x - z)X 
(1 - e* a lx)(3 + e ia l(x - 2z))(E + A) 

( .._ .. j „ (x + z)e- iS ° 



H = -— 



2E 



E\ E + X3 + e ta l(x-2z) 
x(l + l(x - z)e La )e- l& - 



E + A (1 - e ia lx){3 + e la l(x - 2z)) 
K = 



2E 

E V E + X3 + e la l(x-2z) 



2E 



E + X 3 + e la l{x - 2z) 



(33) 




where 

/ z? c n? \" / m /■ 

(34) 

Other coefficients (A, C, D, F, J, L and P) vanish. Let us notice that the modulus square of the above coefficients 
depends only on the phase difference a crossing the interface. Similarly to the case of 2SC/2SC junction this quantity 
we consider is gauge independent . 

The probability current for z < and E > A, 2 A is 

j = *U*< = 2M|(1 - \B\ 2 - \D\ 2 ) - \G\ 2 }. (35) 

The analytical result is not given by short expression thus we rather show the current dependence on the energy in the 
Fig. 1 for generic set of parameters. Let us notice that the current starts at energy E = A and then rises linearly up 
to the point E = 2A where there is a jump. This behaviour is expected because below the gap there is no probability 
current and at the energy twice the gap additional current appears from the singlet excitation in the CFL phase. 
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75 inn 125 150 175 200 225 

FIG. 1. Dependence of the probability current ( p5| ) as a function of energy E [MeV] for two sets of parameters: A = 80 
MeV, A = 60 MeV (black curve) and A = 100 MeV, A = 60 MeV (gray curve). 



For comparison let us show the result for transition and reflection coefficients in the QGP/CFL scattering: 



rp9 . 



T , j for£<|A| 

7 < for^>|A| ^ 

for£<|A| 



sph) fOT £>|A| 
for £: < 2|A| 

2C 



3(S+C) 



for£ > 21 Al 



where T A are transition coefficients for quasiparicles in the CFL basis described by the wavefunctions (|58|) given in 
the Appendix C. Let us note that because of the energy conservation there is always only one hole reflected in the 
QGP phase: 



5 _ 
9 

1 E 



T)d JDS 

green ^blue 



= < 



l(E 2 



9|Ap 

i(E-Q 
9(E+C) 

C _l 4B 

9 E+C 

2 , ! 
9 ^ 

I -L 

9 T 9(E 
1 E-C , 



|C||e|)coB(25)- 
2 [Ecos(28) 



9|Ap 



1 

9 E+C 



9 E+j CO < 2S ) 

|C||£|)cob(25)- 
i 



9 E+C 



«) 

1 E-C 
9 E+C 



imT [Ecos(25) 
i^cos(2£) 



S(|C|-|e|)sin(25)] 
-|C|sin(«5)] 

£?(|C|-|e|)sin(2«y)] 
- |C! sin(5)] 



9 E+C 



iorE < |A| 
for|A| < E < 2|A| 
iovE > 2|A| 

iorE < 2|A| 
for|A| < E < 2|A| 
fori; > 2IAI 



(37) 



where |C| = ^4\A\ 2 -E 2 , |f | = y/\A\ 2 -E 2 . Reflect ion and transition coefficients sum up to unity for any value of 
energy E of incoming particle. The coefficients (37) depends on the phase 5 of the gap parameter. However this single 
phase is a gauge dependent quantity. Thus to get the physical answer we have to avarage over the phase angle. The 
sinus and cosinus avarage to zero giving the final answer in the form: 



ri„„ r i — 



^green ^6 



5 

1 , ME-C) 
9 9(E+C) 
I E-C , 4 E-C 
9 E+C ' 9 E+C 

2 

1 I E-C 
9 ^ 9(E+C) 
1 E-C , 1 E-C 
9 E+C ' 9 E+C 



iorE < |A| 
for|A| < E < 2|A| 
iorE > 2|A| 

iorE < 2|A| 
for|A| < E < 2|A| 
forE > 2IAI 



(38) 



For the case of QGP/2SC scattering the results are the same as given by equations ( |l5| , |l6| ). The only difference 
is that instead of Fermi velocity in the probability current ([l^) we have relativistic expression 2pp (the Fermi 
momentum pp is essentialy equal to quark chemical potential [i in our approximation) . Let us notice that the reflection 
and transition coefficients are exactly the same like in the non-relativistic case. 







VI. CONCLUSIONS 



In this paper we consider the general structure of the Andreev reflection between two superconductors 2SC and 
CFL in the high density QCD. We also give the review of the Andreev reflection between QGP/2SC and QGP/CFL 
phases as well as in condensed matter systems. 

The essence of the AR stems from the peculiar behaviour of particles which hit the interface. If the energy of 
the incoming particle from the conductor side is below the energy gap in superconductor then the hole is reflected 
from the interface. The energy and momentum^] are conserved, however, there is an apparent violation of the charge 
conservation. This last effect stems from the fact that the superconductor serves as an infinite suply of Cooper 
pairs. At the microscopic level charge is conserved, of course. The Andreev reflection can be understood as follows: 
the incoming particle takes another particle from the Fermi sea and creates a Cooper pair which dissolves in the 
condensate. Then the hole is left in the conductor with the appropriate kinematics constaraind by energy and 
momentum conservation. 

This pattern is exactly repeated in QCD superconductors. The differences come from the more complicated structure 
of Cooper pairing of quarks. There are more quantum numbers describing QGP phase as well as the quasiparticle 
structure in QCD superconductors are richer. The treatment of the Andreev reflection processes are similar to the 
case of nonrelativistic physics, however, more attention has to be taken for the gauge invariance of the final results. 
It appears that the probability currents are gauge independent quantities but the transition T and the reflection R 
coefficients require more attention. Our effective hamiltonian is not gauge invariant, nevertheless one can trace 
the gauge behaviour by tracing the dependence of the final result on the phase 8 of the gap parameter. Obviously 
the phase 8 itself is not gauge invariant quantity. In the QGP/2SC case the coefficents T and R are gauge invariant, 
however, this is not true for the QGP/CFL interface (|37|). Let us first notice that one should expect the gauge 
invariant answers for T and R coefficents because the reflected holes can be distinguished by their flavor quantum 
number (not a color) . We suggest that the avaraging of the final result over the phase 8 rebuilts the gauge invariant 
answer. This is reasonable. Indeed to perform the Andreev reflection experiment one has to prepare the ensemble 
of the QGP/CFL interfaces. Then each sample has its own unphysical gap phase. Then performing the experiment 
on the unsamble is equivalent to avarage over the phase 8. Finally in the case of 2SC/CFL interface the physical 
quantities depend on the phase difference a crossing the boundary. This is actualy expected result because, similarly 
to f (1) superconductors, the angle a is a gauge independent variable. 

The Andreev reflection has many interesting physical effects in condensed matter systems |l(],[l(],|l7j . In the high 
density QCD situation is more complicated. The only known place one can expect the color superconductors are 
the cores of protoneutron or Neutron Stars. Any subtle effets of Andreev reflection that happen inside these objects 
far away from our laboratories may be never observed. However there is at least one possibility one can imagine. 
The Andreev reflection affects the transport properties at the interfaces. If the protoneutron stars go through the 
first order phase transitions between QGP or/and 2SC or/and CFL phases during their history cooling then one can 
expect the mixed phase to be present in the core of the stars. In this situation the Andreev reflection influences the 
dynamics of the bubbles grow between different phases. These in turn influence the neutrino emmission. Thus the 
time dependence of neutrino luminosity from supernovae can carry the information of the Andreev processes which 
took place inside the protoneutron stars. These intriguing possibilities require more attention. 

Acknowledgement. M.S was supported by a fellowship from the Foundation for Polish Science. M.S. was also 
supported in part by Polish State Committee for Scientific Reasearch, grant no. 2P 03B 094 19. 

APPENDIX A - Particles, holes and Dirac algebra 

Let us consider the equations of motion for free quarks at chemical potential fj,. 

iip(t, r) = (—id? • V + m7o — m)^(*> r)> 

ift(t,r) = -iVip*(t,r) ■d-<ipi(t,f f )(m'yo - /■*)• (39) 

Remembering that relativistic effects arc suppressed at high density the Dirac field operator ip(t,r) can have a one- 
particle interpretation of annihilating a fermion of given spin at the point (t, r) in space-time. We can also change 



5 The momentum conservation is violated at the 0(1/Ef) level and can be neglected in the first approximation. 
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the operator equation ( p9| ) into the wavefunction equation. The first equation of (|3^) describes the particles and the 
second one the holes. Using the decomposition: 

ip(t, r) = 2J OL r tp r , R {q) exp(iq- r- iet), (40) 

r 

i>Ht,r) =^j*hl L (-q)exp(iq- f- iet), 



one can easily solve the equations (p9|). To finish this let us define the bispinors ip rtR (k) and (fc) through the 
equations: 

(a ■ k + 777,70 - (j,)tp r<R (k) = e(p riR (k) (41) 



where r describes the spin and momentum fc can be in general complex vector. From this reason one has to distinguish 
between the right- and left-handed eigenvectors, which are denoted by the capital letters L, R. For complex k ip\ L '\s 
not hermitian conjugate to tp r ,R- The solutions of the equations take the form: 



/ Jm + e + Zrvv \ 

*V,it(*) = gj (42) 



v4,l(*0 = xJV™ + e + (j,,xt 



<7 • fc 



Xr 



^/m + e + \x 



where Xri xt are spinors and where e = \/k 2 + m 2 — The spinors can be defined in the helicity basis: 

? • fcXT,i = ± fc XU ( 43 ) 

x\ d CT-k = ±k X \ d 

where fc = V^fc 2 ". Let us also define the additional bispinors: 

hl L (k)(a -k-m l0 + fi) = eh\ L {-k) (44) 
(a ■ k - mj + fJ,)h r>R (-k) = eh r>R (-k) 

where e = — e and bispinors are given by formulae: 



/ Jrn — e + uv r \ 
ftr,il(-*) = (45) 



y/m—e+fi 



Xr 



h lA-k) = xtV™- e + fj.,-xt- 



a ■ fc 



y/m — e + \i J 

The bispinors defined above fulfill simple algebraic relations which are useful in the calculations: 



fl.L^^R = K L h s,R = 2\/k 2 + m 2 5rs (46) 



vl,L C l5 fl i'?L h »,R = ( Pj,R C 'Y5h r , R = 2\/k 2 + m 2 | x 1 * = j ^ = | 



The bispinors defined above together with ( fto| ) are solutions of (39) and have simple physical meaning. The 
wavefunction: 



There is another solution with e = + m 2 — /i which is not interesting for our purposes 
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tp(t,f) = <P]B{k) exp (-iet + ik ■ r) (47) 



describes the particle of spin projection up, velocity v = 4, where E = + m 2 and energy e above the Fermi Sea. 
From the other hand the wavefunction: 

i>Ht,r) = h\ L (-k)exp(-iet + ik-f) (48) 

describes the hole of spin projection down, velocity —v, and energy e below the Fermi Sea. 

APPENDIX B - 2SC quasiparticles 



From the hamiltonian ( pL T[ ) and the gap structure (21) for 2SC phase (|23|), using usual commutation relations for 
fcrmi fields, one can derive the equation of motion: 

#red = H« • V + m 7o - fi)^ ed - AC 7s ^ ee „, 

i^green = ■ V + TO 7o - (J<)lpg reen ~ AC 75 l/Ci ( 49 ) 

where we exchange the field operators by the wave functions similarly like in the nonrelativistic case. This can be 
done because the presence of the Fermi sea suppresses the contribution of antiquaries. There is similar set of equations 
for another (u gree n,dred) pair. It is convenient to use the plane wave decomposition: 

Credit, r) =^2a r ifir, R (q)exp(iq- f- iEt), 

r 

^greenit, f) = ]T Prrf, R (q) exp(iq ■ f- iEt) (50) 
r 

where the bispinors tp and h are defined in Appendix A. ot r , (3 r are some constants. The subscripts u, d and s describe 
flavor and color of quarks. Plugging ( |50"| ) into (f49|), using the relations given in Appendix A one can find the wave 
functions [pr| : 

»(M0 S ( & ) = ( expW2)v^(^-CVr,) , ^ 



cM-i6/2)^±f(Dhf L T + D>h*l T ) ) 



ip2Z—iEt 



where S is a phase of the gap parameter, pi.2 = y/JjTi: A) 2 — to 2 , A = \J E 2 — |A| 2 and C,C',D,D' are some 
constants. The first wavefunction describes the particle-like excitation whereas the second one describes the hole-like 
quasiparticle. This is quite similar to the description of the nonrelativistic quasiparticles. Indeed if one compares ( |5l| ) 
with the wavefunctions (10) the similarity is striking. It also occurs that the Andreev reflection in the QGP/2SC case 
is a direct analog of the Andreev reflection at the nonrelativistic conductor/superconductor junction. 

The extension of the wave function (^lj) to the case of the 2SC superconductor with 3 flavors is strightforward. The 
exact expressions are given in Section V. 

APPENDIX C - CFL quasiparticles 

Let us now discuss the CFL quasiparticle waveufunctions for the states u re d, d green , suue- The equations of motions 
that follow from ( |l7| ) take the form: 

^red = • V + TO 70 - /X)^red ~ ^^(^een + KLe), 

dgreen = H« • V + TO 7o - fi)^ reen ~ AC l5 {^*ed + ^blue), 
inline = H« • V + TO 70 - Mine ~ AC 75 (Ced + dgreen), 

*€ld = • « - CU™7o - M) - A*« T ee „ + ^L)^ 75 , 

= -^ 6 S L ' « - ^L(m 7 o - M) - A* (^ + < T ee „)C 75 . (52) 
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To find the quasiparticle wavefunctions for A = const, it is convenient to use the following decompositions: 
Ced(*> ^)=zl a r^, R {q) exp(iq- r- iEt), ip^l d (t, r) = ^2 <4K,l(-Q) exp(i<f ■ f - iEt), 

r r 

^reen(t, = ^ Pr^Aq) exp(iq ■ f - iEt), < een (i, r) = £ fih^ L (-q) exp(iq ■ T- iEt), 

r r 

^biuei^r) = J2^r, R {^ ex P( i 9-r-iEt), ^ lue {t,r) = ^-y*h^ L (-q)exp(iq- r- iEt), (53) 



where ip and h are defined in Appendix A. a r , (3 r and j r are some constants. The subscripts u, d and s describe flavor 
and color of quarks in obvious way. 

Plugging d53| ) into d52| ) , using the bispinor algebraic relations given in Appendix A and assuming constant value of 
the gap parameter A, one obtains the wavefunction describing the quasiparticle excitations of given energy E in the 
CFL phase: 

/ u \ 

I A- r\R \ 



*(t,r) 



/ Wed \ 






ib d 

green 






^red 




H 


green 






V tf e / 







A* rtR 


-h uiT 



"IL 





e iqi ' r + J 



I A* r]R \ 

A* r]R 





-L 



I A' ^T-R » 


A* VjR 



\ 



e il2-r + N 



J 



V 

/ iBc^R \ 
^A_w s 







e il2-r +K 



P 



V 

E+C^IR 





A* <PjR 
' L IL 





"Li 



Mi 



cxp (54) 



where £ = a/P - |A| 2 and C = ^ 2 ~ 4 I A I 2 - 9i,a = a/0" ± £> 2 = ™ 2 and 2>i,2 = V(M±C) 2 -m 2 - ff, J, if, L, AT 
and P are arbitrary constants. Similar expression is obtained for the opposite spin content. Instead of using 
u re d, dg reen , suue basis one can use the CFL-basis: 



A=l VZ 



(A )ia.A 



in which the gap matrix (pll) is diagonal, where 



A , A = 1, 8 are Gell-Mann matrices and A 



Explicitly: 



V> 3 = ^(Ced-^reen) 
^ 8 = 4(Ced + <e en -2^„ e ) 

In the CFL-basis one can write quasiparticle wavefinctions as: 



i>CFL 



particle 
^hole 



4. 



CFL 



v 8 



particle 
^hole 



L 



CFL 



particle 



2H + K ( ^^r 
V2 { ~h\l 



2A 



V3P E +S 



-/i tT 



-/l tT 
"XL 



(55) 

(56) 
(57) 



(58) 
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The excitation v^ CFL describes the singlet quasiparticle with the gap 2 A wheras ^cfL'^cfl describe octet quasi- 
particles with the gap A. The probability currents: 

J A = IPplticle^particle + ^tL^hole (59) 

are conseved and can be used to calculte the transition coefficients in the scattering procesess. 
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